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We design a plasmonic fiber waveguide (PFW) composed of coaxial cylindrical metal-dielectric 
multilayers in nanoscale, and constitute the corresponding dynamical equations describing the modes 
of propagation in the PFW with the Kerr nonlinearity in the dielectric layers. The physics is 
connected to the discrete matrix nonlinear Schrodinger equations, from which the highly confined 
ring-like solitons in scale of subwavelength are found both for the visible light and the near-infrared 
light in the self-defocusing condition. Moreover, the confinement could be further improved when 
increasing the intensity of the input light due to the cylindrical symmetry of the PFW, which means 
both the width and the radius of the ring are reduced. 
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How to control the propagation of the light is a most 
important subject in optics. Using the technology of the 
optical fiber waveguides (OFWs) to pilot the light is a 
big advance towards all-optical signal processing. When 
consider ever-accelerated miniaturization of optical de- 
vices, however, the conventional OFWs seem to be difh- 
cult to fulfill the requirement because of the diffraction 
limitation for the optical components of dielectric pho- 
tonic materials. Recently, it is shown that the limita- 
tion may be overcome in the rapid developing field of 
plasmonics [iLSIl based on properties of the surface plas- 
mon polariton (SPP), a confined mode localized at the 
interface of metal and dielectric material 0, 0] . Subse- 
quently wide attentions have been paid to miscellaneous 
nanostructures with metamaterials involved in pursuit of 
subwavelength confinement of the light [6] . Among these 
progresses, several researchers have predicted the sub- 
wavelength control of the light in various lattices through 
the formation of solitons when the nonlinearity is con- 
sidered. It has been found that the planar nonlinear 
metal-dielectric multilayers (MDM), a stack of alternat- 
ing metal-dielectric nanolayers, can effectively manipu- 
late the propagation of light ''2-|9[ and put a subwave- 
length confinement on it jlOHl5| . For instance, Zhang's 
group theoretically found the subwavelength discrete soli- 
ton in the nanoscaled periodic structures consisting of 
MDM ij] . Ye et al. predicted the stable fundamental 
and vortical plasmonic lattice solitons of subwavelength 
extent in arrays of metallic nanowires embedded in a non- 
linear medium 15 ] . 

Then could we construct a kind of waveguide in anal- 
ogy to the optical fiber waveguide to control the light's 
propagation in two-dimensional (2D) transverse space 
based on the plasmonics rather than the optics? In this 
Letter we design a plasmonic fiber waveguide (PFW) by 
rolling the MDM to a cylindrical shape to guide the light. 
The nonlinearity in the dielectric layers are employed to 
realize the subwavelength confinement in radial direction 
when propagating along the axial direction. Based on 



the advanced nano-technologies in addition to the 

sophisticated fiber fabrications, this kind of PFW is sup- 
posed to be a easily-fabricated structures. Actually a 
spherical hyperlen made of the MDM has been realized 
in experiment [Toj. Besides, in contrast to the lattices 
made of planar MDM, the PFW has the cylindrical sym- 
metry, which produces the centripetal confinement in the 
form of solitons. Our results show that the energy could 
be highly confined in the cross section with size far small 
than the operating wavelength. This kind of confinement 
combines the advantages of fiber's long distance trans- 
portation and minimized transverse space caused by the 
cooperation of nonlinearity and SPP's field enhancement 
effect, which is significant for applications of plasmonics 
in future communication and large-scale integration. 

Our structure of the PFW is composed of alternating 
metal and dielectric coaxial cylindrical layers, as shown 
in Fig. 1(a) and (b), where the coaxial cylindrical metal 
layers are labeled as n — 1, 2, • • • along radial direction. 
The thicknesses of the metal layer and the dielectric layer 
are represented by w and d, respectively. Note that the 
dielectric material at center is a cylinder with the radius 
of d. The dielectric constants of metal and dielectric ma- 
terial are eoCm and coed, respectively, with eo is dielectric 
constant of the vacuum. A set of typical parameters used 
in the Letter are that: the metallic dielectric constant is 
Johnson's data for the metal Ag [i^l, the dielectric con- 
stant for the dielectric material is 10, and their thick- 
nesses are w — lOnm and d — 40nm, respectively. In 
fact, the following theory is still applicable when the pa- 
rameters d and w are functions of n. 

When the light wave with the frequency uj propa- 
gates in the PFW, the SPP modes are excited at each 
interface [2lj|. Therefore the PFW can be considered as 
an plasmonic multiwaveguide system, of which the i-th 
single plasmonic waveguide (SPW) is constituted by the 
corresponding i-th coaxial metal sheathes and the dielec- 
tric environment. The interaction among the propagat- 
ing modes of SPPs in these SPWs would cause the energy 
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FIG. 1: (a) Cross section of the PFW; (b) Schematics of 
the PFW; Only first few metallic layers are draw for demon- 
stration, (c) Dispersion of SPP mode in the 1-th SPW. Ez 
distributions in different SPWs in case of A = 700 nm (d), 
and A = 1550 nm (e), respectively. For clarity, only that in 
SPWs with indexes 1, 2, 6, and 7 are draw. 

diffraction if no other mechanism counteracts it, as shown 
below. We then start by exploring the guiding modes in 
each SPW to gain some insights. By exactly solving the 
Maxwell equations, we can get the dispersion relation and 
field distributions [E^^^^r), i?(")(r)] of the SPP mode in 
n-th SPW. Here we only consider the fundamental TM 
mode of SPP with nonvanishing components , -ff^"^ 

and Er , which are functions of the radial coordinate 
r. The dispersion relation of 1-th SPW is given in Fig. 
1(c), where /3 is the propagation constant of the SPP, 
rib is the refractive index of the environment, and fco is 
the wavenumber in vacuum. The dot-dash line labeled 
by rifj/co is the dispersion of the light in the background. 
From Fig. 1(c) one can see that the propagation constant 
P has a strong dependence on the frequency and the real 
part is always bigger than that of the background, which 
is the characteristics of the SPP. Calculations show that 
/?„ have a small variation for a given frequency, especially 
for the SPWs with high index, which is reasonable as the 
high-index SPWs can be approximated as planar waveg- 
uides with the same width of w. The field distributions 
of the SPP modes in different SPWs is sensitive to the 
frequency. For example, i^i"'' (ri = 1, 2, • • • , 6, 7, • • • ) are 
plotted in Fig. 1(d) for A — 700 nm and in Fig. 1(e) 
for A = 1550 nm. It is clear that the interaction between 
adjacent modes in case of A = 700 nm is small enough to 
take a nearest-neighbor approximation while it is not in 



case of A = 1550 nm. 

For the whole PFW, we can express the total elec- 
tric field E{r, z) and magnetic field i?(r, z) in the 
superposition of the modes in all SPWs: Er — 

E„a„(^)£^r 'e-^"*, E, = Y.nan{z)—Er e-^^' , and 

— J2n'^niz)H''^'' e~^'^* , where a„ is the mode am- 
plitude of the n-th SPW, e and e^"^ are the dielectric 
constant of the PFW and of n-th SPW, respectively. For 
simplicity all these modes of SPWs have been normalized 

by their energy flow In = ^ Jg Re (-E^"^* x ff^")) • e^dS 
here. Thus the equation describing the dynamics of the 
mode amplitudes a„ can be determined by the general- 
ized Lorentz reciprocity theorem p^. which can be writ- 
ten as: 

«4cA-f (BC-K)A-r|ApA = 0, (1) 

dz 

where A is a vector given by the elements a„ (n = 
1, 2, • • • ) and | Ap represents a diagonal matrix with cor- 
responding elements |a„p. Matrix B is a diagonal matrix 
with the elements given by the propagation constants /3„. 
C, K and F are matrices, whose elements are defined re- 
spectively by 

C„,„, ^ \ f {e^"'^H^;^ + Ei-^H^^) dS, (2) 

if„,„, = ^^(e(")-e)i^„,„,d5, (3) 

r„,. ^ Z!2i|^i^„ „,d5, (4) 

("') 

where K^,, = Ei'^'^Ei'^^ - t—E'f'^E'f^ , a„ = 

|£;(«)|2 + |l^^(")|2^ 

In the above derivation for Eq. ([ij , the change of rel- 
ative dielectric constant td due to the introduced non- 
linearity in the dielectric layers is ^e^|i?p with the Kerr 
coefficient ^. Also, we have neglected the nonlinear in- 
teraction among the SPWs, i.e. items including |a„'pa,i 
{n 7^ n') are neglected. If the distance d is big enough 
so that the overlapping of adjacent SPWs' fields can be 
ignored, the matrix C can also be approximated as a di- 
agonal matrix, as shown in Ref. [l^ . But for the case of 
strong interaction, as the case shown in Fig. 1(e), the ap- 
proximation is not acceptable. Therefore the Eq. ([I} is 
suitable to more general cases, which can also be written 
as a more simplified form 

ii-A + TA + GIApA^O, (5) 

dz 

where T = C^l^C - C^K and G = -C-^r. This is 
discrete matrix Schrodinger equations. For the param- 
eters we adopted, we find that matrix T and G can be 
approximated as a triple diagonal matrix and a diagonal 




FIG. 2: (color online) (a) the profile of the amplitude of the 
longitudinal field component of the soliton with the wave- 
length of A = 700 nm in the PFW. (b) the distribution of 
along the radial direction, inset is an amplification of one 
part for clarity, (c) The nonlinear (soliton) and (d) the linear 
propagation over 40 ^m distance in lossless PFW. (e) The 
propagation of the soliton in the same PFW when the loss in 
metal is taken into account. 



one, respectively. Neglected elements are at least one or- 
der of magnitude less than others. Then it is reduced to 
the nonsymmetric discrete nonlinear Schrodinger equa- 
tions with variable nonlinear coefBcients. In the equa- 
tions, the matrix T is determined by the interaction be- 
tween the SPP mode in each SPW with with its coun- 
terpart in the inner nearest SPW and outer one. Due to 
the cylindric symmetry of the PFW, T is nonsymmetric 
matrix, which leads to the different diffraction towards 
the center or outside. The matrix T governs the nonlin- 
earity, which will balance the diffraction when solitons 
forms. The elements \Gn\ [n = 1, 2, • • • ), proved to de- 
crease with the n index, give a centripetally increasing 
nonlinearity. Both these characteristics contribute to the 
highly confined solitons in the PFW. 

We sought the soliton solutions in the form of a„ {z) = 
•\/7ou„e'''^, where the amplitudes u„ and p are both in- 
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FIG. 3: (a) The propagation of soliton along 40 ^m in the 
lossless PFW when the wavelength A is 1550nm. (b) Linear 
propagation along 4 /im in the same PFW. (c) The propaga- 
tion of the soliton along 4 /^m when the loss is considered. 



dependent of z. We also define the intensity as / = 
in the solution. We put our emphasis on the 
unstaggered solitons [13] formed in the self-defocusing me- 
dia, with the Kerr coefficient set as — 1 x lO'^^m^ /V^ and 
Iq about 5 X lO^^VF. The number of layers is chosen big 
enough to ensure the fields vanishing at the boundary. 

Fig. 2(a) plots the profile of the amplitude of the lon- 
gitudinal field component of the soliton in the PFW 
with / = 0.05/o and Fig. 2(b), with an amplification 
of one part, presents the distribution of along the 
radial direction in the case of A = 700 nm. As E^ re- 
flects the intensity of the light, we can conclude from 
Fig. 2(a) that the energy mainly concentrates inside the 
region with the radius of about 400 nm. It is noteworthy 
that the maximal intensity is at the circle with radius 
about 250 nm other than at the center, which could be 
seen it more clearly from Fig. 2(b). That is why we 
call it as the ring-like solitons. It is the nonlinearity that 
confines the energy in a region with the size of subwave- 
length, where the effect of field enhancement provided by 
the SPPs is prominent, which in turn makes the nonlin- 
earity more achievable. The cylindric symmetry enable 
energy to concentrate towards the axis, which further fa- 
ciliates the confinement of energy. The linear case where 
the diffraction dominates are also plotted in Fig. 2(d) for 
comparison. Moreover, Fig. 2(b) also shows that metal 
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FIG. 4: (a-c) The profiles of the sohtons with the intensity 
of 20/o, lOO/o, I8O/0 respectively. The values in each figure 
have been normalized to their maximum, (d) Radius of the 
soliton's peak changing with the intensity. 



layers hold a large portion of the energy. For this reason 
in addition to a relative big value of Im[/3„] for visible 
light, the soliton experiences a high loss in the PFW if 
the loss in the metal is considered, as shown in Fig. 2(e). 
A gain medium is consequently necessary for the PFW 
working on the visible lights. 

In the case of A = 1550 nm, the ring-like solitons are 
also found in the same PFW. Its propagation along 40/im 
in the lossless situation with / = 8O/0 is plotted in Fig. 
2(a), while the lossless linear propagation is in shown 
in Fig. 2(b). Compared with the case of A = 700 nm, 
the diffraction is outstanding, because of the remarkable 
overlapping between the adjacent SPWs' fields and thus 
a strong interaction among the SPPs in different SPWs, 
as we discussed before. And consequently the highly con- 
fined soliton is achieved with a relative strong intensity. 
However, as the portion of the field residing inside the 
metal layers decreasing correspondingly, the transmission 
distance is enlarged to some extent in the lossy PFW, 
which is clearly seen by comparing Fig. 3(c) with Fig. 
2(e). 

For the lattices solitons, increasing both the intensity 
and the nonlinear coefficient will enhance its nonlinear ef- 
fect and therefore make the energy more confined. In the 
PFW, this characteristics is particularly obvious. The 
profiles of the solitons under three values of the intensity 
/ = 20/0,100/0,180/0 are drawn in Fig. 4(a)-(c). The 
energy of the soliton is confined in a ring of less and less 
width with the increasing of intensity. Furthermore, be- 
cause of the nonsymmetry of the matrix T, increasing 
the intensity also makes the energy moving centripetally. 
It is an exciting result that the highly confined 2D soli- 
ton in Fig. 4(c) has a transverse size far small than the 



wavelength. The peaks of these solitons with more values 
of intensities are plotted in Fig. 4(d). Calculation shows 
that when the intensity reaches a certain value, the peak 
stops moving. The reason is that in this situation the 
nonsymmetry of the diffraction is trivial for the consid- 
erable nonlinearity. Besides, as the matrixes T and G 
are determined by the field distribution of SPWs, the 
soliton shape could be modified by adjusting the radius 
and width of each metal layer. 

In conclusion, we designed a PFW composed of the 
MDM and predicted theoretically the form of highly con- 
fined subwavelength solitons in the PFW with Kerr non- 
linearity both for the visible light and the near-infrared 
light. The equations describing the solitons have been 
obtained in an compact matrix form. This kind of 
PFW is expected to take the role of OFW in future in 
the plasmonics-based communications and other nano- 
photonic applications, such as lithography, beam shap- 
ing, et ai. 
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